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The concept of Carnot efficiency is a founding principle of macroscopic 
thermodynamics. It allows to introduce entropy as a state function and to define 
the Kelvin temperature scale. It states that, for a system operating between two 
reservoirs at temperatures and T c , the efficiency fj = W/Qh , being the ratio of 
output work W over input heat Qh , is bounded by the Carnot efficiency fj < rjc with 
r/c = 1 — T c /Th. Right from the start, the question was raised about the efficiency 
of small scale machines. Maxwell introduced a small-scale demon which was deemed 
to rectify thermal fluctuations. A clarifying rebuttal was given by Smoluchowski, who 
proposed a mechanical implementation of the Maxwell demon, the so-called ratchet 
and pawl. He stressed that this small-scale device would eventually thermalize, after 
which any rectification would stop. Szilard introduced an information driven engine 
which seemed to be able to extract work from a single reservoir, in violation with the 
Carnot prediction. Since he believed that this could not be true, he concluded that 
there had to be a thermodynamic cost associated to the information gathering [1], 
The Szilard engine gave rise to a prolonged scientific discussion about the source of 
dissipation [2], It is not the measurement or computational process, but the resetting 
process of the memory device that seems to be the dissipative step. The Smoluchowski 
engine was revisited by Feynman, who showed by an explicit model calculation that 
the efficiency of the ratchet and pawl in contact with two reservoirs is indeed bounded 
by the Carnot efficiency. A subtle error in his analysis was elucidated by Parrondo 
and Sekimoto [3, 4], indicating that the efficiency was strictly below Carnot efficiency 
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in this model. This was confirmed by a more detailed analysis on simplified models 
[5, 6], prompting the question whether Carnot efficiency could at all be reached in such 
small devices [7]. These questions have led to major interest in information to work 
conversion, both theoretically [8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18] and experimentally 
[19, 20, 21, 22, 23, 24, 25]. 

Thermodynamic efficiency can also be defined for other types of engines, notably 
for work to work and the above mentioned information to work transformations. In 
the case of a transformation of input work Wi to output work W a , the thermodynamic 
efficiency is defined as fj = W Q /Wi. The analogue of Carnot efficiency is reached for a 
reversible operation leading to W Q = Wi (since no heat is dissipated). Hence the second 
law stipulates fj < r] r , with the reversible efficiency is r/ r = 1. For information to work 
transformation, we note that information about the system can be used to extract work. 
The reversible limit has been known since the work of Szilard, namely k^T In 2 of work 
can be extracted per bit of information (in an environment operating at temperature 
T ). One bit corresponds to a Shannon information of / = In 2. More generally, in a 
transformation of a Shannon information amount I into an amount W of work, the 
efficiency f) = W/(ksTI ) is upper bounded by the reversible result r/ r = 1. 

Over the past two decades, one has been able to reformulate thermodynamics 
to describe fluctuating small-scale systems. The most notable of the results is the 
fluctuation theorem, stating that the probability for a stochastic entropy production 
A iS is exponentially more probable than that of a corresponding decrease —AjS in an 
’’inverse” experiment: 

_j0_ 

P(-\s) 

The tilde refers to the “time- inverse” experiment. The above symmetry relation for the 
probability implies the following integral fluctuation theorem: 

(exp (—fceAjs)) = 1, (2) 

which in term implies the “second law”: 

<Ais) > 0, (3) 

The detailed and integral fluctuation theorems (1) and (2) have been derived in many 
different settings, see for example the Jarzynski [26] and Crooks [27] equalities, and 
stochastic thermodynamics [28, 29]. When considering a finite time experiment, the 
interpretation of (1) is somewhat delicate, and there may be restrictions, for example 
on the initial and final states of the experiment, cf. [30] for a recent discussion. The 
fluctuation theorem however appears to have a wide ranging validity, comparable to 
the second law, when considering the asymptotic long time limit, in which case (1) 
reduces to a statement about the large deviation properties of the entropy production. 
The theorem was actually first derived in this context [31, 32, 33, 34], Surprisingly, 
the impact of these new insights on the efficiency, and in particular on its stochastic 
properties, has only been considered very recently [35, 36, 37, 38, 39, 40, 41, 42, 43], 


exp(fceAjs). 


(1) 
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and only in work to work and heat to work converters. When running a small-scale 
engine for a finite time, the corresponding cumulated work output and heat uptake, 
w and qh, or work output w a and work input w*, are stochastic quantities, and hence 
so is the corresponding stochastic efficiency rj = w/qh or rj = w Q /wi. Starting with 
the detailed fluctuation theorem ( 1 ), it was pointed out that one can make universal 
statements about the stochastic efficiency 77, as it approaches the macroscopic efficiency. 
More precisely, for large but finite times, values of the stochastic efficiency rj different 
from the macroscopic efficiency fj are exponentially unlikely, as described by the large 
deviation function J[rj) = — linp^oo l/t ]nP t (rj), J{rj) > J{fj) = 0. For the case of 
a thermodynamic machine driven by a time-symmetric protocol, it was shown that 
the reversible efficiency is the least likely, i.e., the large deviation function J(rj) has a 
maximum at 77 = q r . The macroscopic efficiency is reproduced as the value carrying 
all probability to dominant order J(fj) = 0. For time-asymmetric protocols, the large 
deviation curves of the stochastic efficiency for the forward and backward experiment 
cross at the reversible efficiency [36, 37]. 

In this paper, we provide an explicit and comprehensive illustration of the stochastic 
efficiency in five engines. Providing possibly the simplest steady state models for work 
to work and heat to work transformation, we consider a Brownian particle subject to 
competing forces [35] and effusion between two compartments [38]. Even though these 
cases have been discussed in some detail in the literature, they are introduced briefly for 
completeness and for comparison with the other models. The third model is a thermal 
engine based on a quantum dot. It is of interest because its stochastic thermodynamic 
properties have been discussed and measured, and because it is in principle richer than 
the effusion model, which can be recovered in an appropriate mathematical limit. In the 
last two problems, we evaluate stochastic efficiency in novel settings. We show that the 
universal features of stochastic efficiency are valid for heat to momentum transformation. 
This is illustrated on an effusion model with momentum transfer. Finally, we discuss 
the stochastic efficiency for information to work transformation. It is indeed possible to 
reformulate the fluctuation theorem when dealing with an information processing set¬ 
up, such as the one introduced by Szilard. In particular one can introduce the stochastic 
efficiency 77 = w/i for a machine transforming a stochastic amount of input (Shannon) 
information i into work w. We show that its large deviation function again displays the 
same general features. In particular, the reversible efficiency q r = 1 is exponentially 
less likely than any other efficiency in the asymptotic time limit (for time-symmetric 
protocols). We illustrate these features on the Mandal-Jarzynski model [8, 10]. 

1. Brownian Engine 

Consider an overdamped Brownian particle on a plane, subject to two external forces, 
a loading force iq, and a driving force F 2 [35], cf. Fig. la. The (larger) driving force 
F 2 pushes the particle against the loading force F L . The stochastic efficiency of such a 
device as a work-to-work converter was discussed in [35, 39]. The mathematics are very 
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(a) Schematic representation of a Brow- 
nian work-to-work converter 



(b) Colour-coded macroscopic efficiency 
of a Brownian engine in function of 
F 2 with Fi = (1,0). The green dot 
corresponds to F 2 = (—3/2,1) and fj = 
0.29. The stochastic efficiency for this 
case is represented in Fig. 2a. 


Figure 1: Characteristics of the Brownian engine 


simple. Considering for simplicity a two dimensional set-up, the displacement x of the 
Brownian particle during a time t is characterised by a bi-Gaussian distribution. Under 
influence of the resulting force F = F\ + F 2 the average displacement is (x) = / iFt , /i 
being the mobility. The variance is isotropic and uncorrelated in orthogonal directions, 
(SxSx) = 2 Dil, where D is the diffusion coefficient and 1 the unit matrix. 

We first turn to the macroscopic efficiency, which is very easy to evaluate, see also 
Fig. lb for a colour-coded illustration: 

n - W - FSF = FiF (4) 

<u>2> F 2 -(x) F 2 -F' 

where W\ = F\ • x and w 2 = F 2 ■ x are the stochastic amounts of work delivered by the 
loading and driving force respectively. The engine regime, i.e., the regime where the 
driving force delivers a positive amount of work to the loading force, is determined by: 


< cos 9 < 




( 5 ) 


with 6 is the angle between F\ and F 2 . In combination with (4), it is clear from (5) 
that, in the engine regime, the macroscopic efficiency is bounded by fj < rj r — 1 . The 
reversible efficiency r) r = 1 can only be reached in the limit F 2 —» —F\ with F 2 || —F±. 
This can also be seen in Fig. lb. 

We next investigate the stochastic efficiency rj = —w\jw 2 = —F- v -xjF 2 -x. Being the 
ratio of two correlated Gaussian variables, its probability distribution can be evaluated 
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(a) Probability distribution P t (r 7 ) of the 
efficiency 77 for the Brownian engine. 


5 



(b) Approach of — \nP t (r])/t for in¬ 
creasing values of t/to (blue, red and 
green curve) to the large deviation func¬ 
tion J(? 7 ) (black curve). The yellow 
curve is obtained by extrapolation from 
the finite time results. 

= (1,0), F 2 = (-3/2,1) 


Figure 2: Efficiency fluctuations of a Brownian engine with Fi 
and fj = 0.29. 


analytically, see also [39]: 

F x x F 2 


«(») = 


e *o 


(F 1 + V F 2 )H 


1 + V 7 r 5 '( r ?)Erf(\/ 5 r (??))e 9( ' ?) ) , 


( 6 ) 


with 


g{v) = 


t (1 — 77 ) x F 2 
F 2 (Fx+r^Fa ) 2 


to 


(7) 


determines 


and Erf(x) is the error function. The characteristic time to = 2 D/ y/j, 
the boundary between diffusion dominated ( t -C to) and drift dominated (t to) 
dynamics. We note in passing that it is easy to show from equations ( 6 ) and (7) that 
pm > 0 and F/(r^ r ) < 0, implying that there exists at least one maximum in the 
interval rj £ [ 0 , 1 ] 

Universal features of the efficiency fluctuations are revealed when studying the 
asymptotic time behavior via the large deviation function of r] [35]: 

1 


J(rj) = — lim - In Ptijf) 

t —>00 t 


/' 


(Fi + rjF 2 ) ■ (Fi + F 2 ) 


1 2 


4F> 


(Fx+r ?^) 2 


( 8 ) 






















Stochastic Efficiency: Five Case Studies 


6 




Uh ! 'i)i fic , T c 

(a) Schematic representation of the 
effusion model. 


0 



0 
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(b) Macroscopic efficiency of the effu¬ 
sion engine in terms of UhlksTh and 
Hc/kBTh, with r/c = 1/2. The green 
dot corresponds to Uh/ksTh = — 1 and 
Hc/kBTh = —3/4. The stochastic ef¬ 
ficiency for this case is represented in 
Fig. 4b. 


Figure 3: Characteristics of the effusion model 


This function has a minimum at the macroscopic efficiency J(fj ) = 0, but also a 
maximum at the reversible efficiency r/ r = 1. with equal asymptotes in the limits 
ri —» ±oo. To illustrate the approach to the large deviation regime, — ]n(P t (rj))/t is 
plotted in Fig. 2b for t/to = 2,5, and 10, together with the limiting expression (8). 
We also include the result of an extrapolation ansatz [38], described in more detail in 
Appendix A. The extrapolation, based on the t/t 0 = 2, 5, and 10 curves, is in surprisingly 
good agreement with the exact asymptotic expression. Although these finite time results 
do not, in this particular instance, exhibit a maximum close to the reversible efficiency 
Tj r = 1, it does show up by extrapolation. 

2. Effusion Engine 

The effusion engine [38] consists of two reservoirs, exchanging heat and particles by 
effusion via one or more small holes in the separating wall, cf. Fig. 3a. The reservoirs are 
supposed to be infinitely large and at equilibrium. The holes are smaller than the mean 
free path so that the equilibrium state is not disturbed by the effusion process. Under 
proper working conditions, a net flux of particles moves from say the left compartment, 
at high temperature T/ and low chemical potential /.(/,, to the right compartment at 
lower temperature T c but higher chemical potential // c . 

When a particle moves from the hot to the cold reservoir, it delivers an amount of 
work Wq = fi c — jin = A//, while extracting an amount of heat <jq = u 0 — ///,, from the 
hot reservoir, where uo is the kinetic energy of the transfered particle. Therefore, after 
a net transfer of n particles, the total amount of delivered work and extracted heat, w 













Stochastic Efficiency: Five Case Studies 
and q are given by: 


7 


w = nAfj,, 


( 9 ) 


q — u — niJL h , 


( 10 ) 


with u the net energy transfers. 

Particles with kinetic energy E transfer from the hot to the cold reservoir at rate 

[44]: 

T^c(E) = ( 11 ) 

to ( k B T h ) 

and from the cold to the hot reservoir at rate 

1 E 


T c ^ h (E) = 


_ E | Me _ ffh 

k B T c ' k B T c k B T h 


( 12 ) 


*o (k B T h f 

where t 0 = \J‘It: rn / (a 2 p/ /.: R E ,) is the average time between particle crossings from the 
hot to the cold reservoir, a is the surface area of the effusion hole and m is the mass of 
the particles of the gas. The macroscopic efficiency can be easily obtained: 

(w) 


V = 


(q.) 


A /j, (n) 


(u) - ia h ( n ) 

A fi (ik B T h ) 2 e*tfc - ( k B T c ) 2 e 




{ksThY^ksTh — iJ,h)e k B T h — ( ksT c ) 2 (2 ksT c — /pj e 


Me 


(13) 


This is plotted for Carnot efficiency r]c = 1/2 in Fig. 3b in terms of Hc/kuTh and 
[i h jk B T h . The engine boundaries are given by [i h < // c < fi h T c /T h - 2k R T c hi (T c /T h ). 
Furthermore, the macroscopic efficiency is bounded by the Carnot efficiency, and this 
boundary is only reached in the limit p c , Hh —> —oo, which, for an ideal gas, corresponds 
to zero density. 

We next investigate the stochastic efficiency r] = w/q. Since it is not possible to 
obtain the analytic expression of the probability distribution Pt{rj), we present results 
from a numerical simulation of the Markov process with the prescribed rates (11) 
and (12), cf. Fig. 4a. For the parameter values under consideration, one clearly sees 
a minimum in the probability distribution developing in the vicinity of the Carnot 
efficiency rjc , even at the rather short times represented here. The behavior around 
77 = 0, and in particular the other minimum around i] = 0, can be explained by the low 
number of particle crossings for short times [38]. 

Turning finally to the asymptotic time behaviour, we note that the large deviation 
function J{rj) of the stochastic efficiency can be obtained from the joint cumulant 
generating function of work and heat cp(A,cv), cf. [37]. The latter is explicitly known for 
effusion [44] : 


ld(A,uj) = lim-ln(e xw Ul( ^) 

’ ’ t->oo t x ' 
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(a) Probability distribution P t ( 77 ) of the (b) Approach of —\nP t (rj)/t for in¬ 
efficiency r] for the effusion engine. creasing values of t/to (blue, red and 


green curve) to the large deviation func¬ 
tion J(rj) (black curve). The yellow 
curve is obtained by extrapolating the 
finite time results. 


Figure 4: Efficiency fluctuations of the effusion engine, with rjc = 1/2, ///, = —/crT/,, 
and ji c = —?>/ AksTh. The macroscopic efficiency is given by fj = 0.07. 



(14) 


The iarge deviation function of the stochastic efficiency is then given by: 
J{rf) = — min</?(A, Xr/). 

A 


(15) 


The contraction can be done numerically, and the comparison with finite-time 
simulations is represented in Fig. 4b. Note that the extrapolation again works quite 
well, except in the vicinity of rj = 0. 

3. Quantum Dot 

The quantum dot model, schematically represented in Fig. 5a, has been investigated 
in detail in the context of stochastic thermodynamics [45, 46, 47, 48]. Two electron 
reservoirs are brought in contact with each other via one or multiple quantum dots. In 
order to investigate its stochastic efficiency, we focus on the case of two quantum dots, 
each with one “active” energy level, E\ and E 2 {E\ < E 2 ). Occupancy of a quantum 
dot by multiple electrons is forbidden because of Coulomb repulsion. For mathematical 
simplicity, we also set all coupling constants between dot and reservoirs equal to T. 

The operation of each quantum dot as a thermal engine is similar to that of the 
effusion engine: a net motion of electrons from a reservoir with low chemical potential 






























Stochastic Efficiency: Five Case Studies 


-C 

=L 

1- 


u 

=L 

u 

1- 


E 


/ 

l-^c 



Mh 



(a) Schematic represetation of the 
quantum dot. 
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(b) Macroscopic efficiency of the quan¬ 
tum dot, with ye = 4/5, E x = ksTh 
and E 2 = 2 


Figure 5: Characteristics of the quantum dot 


to one with higher chemical potential is induced by a driving temperature gradient. For 
every particle transferring through the quantum dot with the lower energy level E x , the 
heat taken from the hot reservoir is given by Sq x = E\ — fih and the delivered amount of 
work is wq = n c — Hh = A/r. Furthermore, the rate of transfer between the hot reservoir 
and the quantum dot are given by: 




K 


Ye T h 


1 + e" 


E \-Ph ’ 


(16) 




r 

E 1 ~^h ’ 

1 1 ■ 


(17) 


The rate of exchange between the cold reservoir and the quantum dot is obtained by 
replacing with Th and Hh by T c and /i c , respectively. Analogous expressions hold for 
transfer through the other quantum dot, with E\ replaced by E 2 . The total amount 
of delivered work w and consumed heat q after a net transfer of n\ particles through 
the quantum dot with energy level Ei and n 2 particles through the quantum dot with 
energy level Eo are then given by: 


w = An (ni + n 2 ), 

(18) 

q = Sq x ni + Sq 2 n 2 . 

(19) 


The macroscopic efficiency can now be written as: 

(w) 


V 


(q) 


A/i((ni) + (n 2 )) 


(ni) (E x - n c ) + (n 2 ) (E 2 - /r c ) ’ 


( 20 ) 


( 21 ) 
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(a) Probability distribution P t (■ r )) of the 
efficiency y for the quantum dot engine. 
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(b) Approach of —In P t {r])/t for in¬ 
creasing values of t/to (blue, red and 
green curve) to the large deviation func¬ 
tion J(y) (black curve). The yellow 
curve is obtained by extrapolation from 
the finite time results. 


Figure 6: Efficiency fluctuations of the quantum dot, with r]c = 1/2, E\ = A’bX/,, 
E 2 = 10 ksTh, i-i h = —ksTh, Pc — 0 and to = 10 5 /r. The macroscopic efficiency is given 
by fj = 0.09. 


where rq and n 2 are the (stochastic) net amount of particles transferred through 
respectively the quantum dot with the lower and higher energy level. Their average 
value can be calculated from equations (16) and (17): 


fai) 


T 

2 


■El-Me £ 1 ~M)i. 

e t c - e T h 



( 

le T h 



( 22 ) 


and an analogous expression for (n 2 ) with E\ replaced by E 2 . The macroscopic efficiency 
of the engine is plotted in Fig. 5b in function of Ph/ksTh and iiJUbT),. These results 
are comparable with the results of the effusion engine. Again, Carnot efficiency is only 
reached in the limit of zero density. 

For the study of the efficiency fluctuations at finite time, we again first turn to 
numerical simulations. The probability distribution Pt(rj), rj = w/q, is obtained by 
sampling the net fluxes rq and n 2 using the rates specified in equations (16) and (17). 
A typical result is shown in Fig. 6a. The results appear to be rather noisy, which is 
due to the fact that the delivered amounts of work and heat are discrete variables. 
Nevertheless, the minimum at Carnot efficiency is very striking, even at these short 
times. The convergence to the macroscopic efficiency on the other hand is rather slow, 
which is, in this particular case, due to the small value of the large deviation function. 

Due to the fact that this engine is driven by two independent, tight-coupled 
operating channels (namely the two quantum dots), with heat consumption and 
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(a) Schematic representation of effusion 
with momentum 



(b) Colour-coded representation of the 
macroscopic efficiency for effusion with 
momentum transfer for 77 c = 1 / 2 . 

The green dot corresponds to Vh = 
5 yJk B T h /m, V c = 5.5 ^/k B T h /m and 
fj = 0.57. The stochastic efficiency for 
the latter case is represented in Fig. 8 a 


Figure 7: Characteristics of the effusion engine with transversal momentum. 


efficiencies per particle 5qi and rji = w 0 /Sqi respectively, i — 1 , 2 , the large deviation 
function for the efficiency J{rj) can be written in terms of the event large deviation 
functions <^(n), i — 1, 2 (see Appendix B): 

J(n)= min (* (H) + <h (yy ) ) ' (23) 

Here, a = (77 — ^ 72 )/(^ 7 i — V 2 ), and 4>i(n) and 02 (n) are the large deviation functions of 
the net number of transferred particles through each of the channels. As these large 
deviation functions are known (see Appendix C), equation (23) can be used to estimate 
the large deviation function of r/. The comparison with finite time numerical simulations 
is shown in Fig. 6 b. Note that the extrapolation from finite time results does not work 
uniformly well due to the aforementioned discreteness of work and heat variables. 

4. Effusion with momentum transfer 

The effusion model, discussed in section 2, has also been studied in the presence of 
momentum exchange between the reservoirs [49], see Fig. 7a. In this set-up, the gases 
move with overall average speed 14 and V c parallel to the separating wall containing the 
effusion hole. For the purpose of illustration, we assume equal densities and chemical 
potentials in both reservoirs, and consider a thermal engine, with temperatures 20 and 
T c in the respective compartments, driving momentum exchange. 

Note that there are 3 stochastic fluxes: the transfer of particles, of energy and of 
transversal momentum. This leaves some freedom in the definition of the efficiency. 
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Here, we choose to define it as: 

A N n + A Px p x 
V ~ A v u 


with 





A 

A 


Px 


Vn 

T h 

1 



Vc 

T c ' 

1 

Th' 


(24) 


(25) 

(26) 
(27) 


the affinities of the particle, momentum and energy transport, and n, p x and u the 
amounts of particle, momentum and energy transport. The macroscopic efficiency then 
reads: 

- _ An (n) + Ap x fp^_ 

^ Au ( u) 

with (n), (p x ) and (it), the macroscopic particle, energy and momentum fluxes. The 
total entropy production is given by [49]: 



A ,S = A N (n) + (p x ) + A v (u) > 0. 


(29) 


Hence the second law of thermodynamics dictates that macroscopic efficiency in the 
engine regime is smaller than 1. This is illustrated for the parameter values considered 
in Fig. 7b. 

Following the discussion of the previous examples, we now turn to the efficiency 
fluctuations. For the finite-time probability distribution we again rely on 

numerical simulations. The rates of particle transfer are given by: 


T h ^ c (E,p x ) = 


1 


V2 


to (k B T h ) ity/m 


exp 


m 




E-^-\ 
2m ) 


1/2 


2 (E-t 


m 


+ 


- + v h 

m 



(30) 


T c ^ h (E,p x ) = 


V2 


( k B T c f 2 (k B T h ) 1/2 ny/m 

2 



1/2 


exp 


m 


2 ksT 



(31) 


with t 0 = v / W(fc^V), the average time between particle crossings from reservoir 
h to reservoir c and p the particle density inside the reservoirs. The results of the 
simulations are shown in Fig. 8a. The minimum at the reversible limit is clearly visible. 
Note also that, in contrast to effusion without momentum transfer, the probability 
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(a) Probability distribution P t (r 7 ) of the 
efficiency 77 for effusion with momentum 
transfer. 
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(b) Approach of — \nP t (r])/t for in¬ 
creasing values of t/to (blue, red and 
green curve) to the large deviation func¬ 
tion J(? 7 ) (black curve). The yellow 
curve is obtained by extrapolation of 
the finite time results. 


Figure 8 : Efficiency fluctuations of effusion with momentum, for rjc = 1/2, 14 = 
5 \Zk B T h / m and V c = 5.5 xJkuTfJrn. The macroscopic efficiency is given by r/ = 0.57 


distribution behaves smoothly around 77 = 0. The explanation is that, due to the 
possibility of momentum transport alone, small efficiencies are possible without net 
particle transport. 

To evaluate the large deviation function of the efficiency, we first note that: 

-WA1V— \qAQ\^ __ ^—\ W (A N AN+A Px Ap x )—\ q A u AU\^ (32) 

and therefore the cumulant generating function of the produced work and heat 
Ho(Xw,^q) can be written in terms of the cumulant generating function of the 
transferred momentum, energy and particle numbers /p (Xu, X N , X Px ): 

fi 0 (X w , X q) = fii(A N X w , A Px X w , AjjXq) 


with: 


= <r(^-) 

\ 27 rm J 


1/2 


pVi 

1/2 




1/2 / 1 _ Gh ( AnXw , AuXq ) 

(1 + ksThAuXQ ) 2 

jU /2 j ^ _ G c {A N X W , A px X W , AjjXq) 


(1 — kBT c AuXQy 


Gh (Aat, X Px , Xu) — exp I —A# — 


mVyXu - k B mT h Xp x + 2mV h X 


x Px 


2(1 + k B ThXu ) 


(33) 


(34) 


G c (Atv, X Px , Xu) — exp ^Atv + 


mVyXu + kBmT c Xp x + 2mV c X 


'Px 


2(1 — k B T c Xu) 


(35) 


and 
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(a) Schematic representation of the 
Mandal-Jarzynski model with 3 energy 
states. 



(b) Alternative physical implementa¬ 
tion of the Mandal-Jarzynski model 


Figure 9: Possible Set-ups for the Mandal-Jarzynski model 


As was discussed for the effusion model, the large deviation can then be found by 
numerically contracting the cumulant generating function: 

J(rj) = - min/i 0 (A, Xij) = — mmfii(A N X, A Px \ AjjijX). (36) 

A A 

The comparison with simulations is shown in Fig. 8b. The extrapolation ansatz seems 
to work quite well apart from an overshoot at the reversible efficiency. 

5. Mandal-Jarzynski model 

Recently, there has been considerable interest in the stochastic information-to-work 
conversion. To study this issue in the context of stochastic efficiency, we focus on one 
of the simplest models, namely the Mandal-Jarzynski engine [10], cf. Fig. 9a. 

A particle, in contact with a heat bath at constant temperature T, can make 
transitions between a number of energy levels, which are multiples of SE. Apart from the 
thermal dissipative transitions due to the bath, the particle undergoes transitions that 
are driven by the entries of a linear tape. This interaction corresponds to a form of input 
work. The energy levels of the systems are numbered as shown in Fig. 9a. The input 
tape consists of a sequence of entry values, each referring uniquely to one of the possible 
energies, which are presented subsequently to the system. During such an interaction, 
the particle is moved to the energy level with number equal to the initial value on the 
tape. In this way, work is delivered. After each entry of the tape, the system is left to 
thermalize, and the new energy state of the particle is written into the tape, which we 
call the final or exit value. Note that this will alter the Shannon entropy of the tape. 
The tape moves on one step and the process is repeated. We mention a closely related 
implementation of the Mandal-Jarzynski model, shown in Fig. 9b. A Brownian particle 
in contact with a thermal bath at temperature T, can make jumps on a discrete lattice 
with a wall at position zero and a force gradient in the negative direction, inducing an 
energy difference between two neighbouring sites equal to SE. After the particle has 
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thermalized, its position is measured, and the wall is instantaneously brought to the 
measured position. Next the wall is adiabatically moved back to the original position 
with delivery of work, after which the process is repeated. While we studied both 
versions, we focus in the following on the simplest situation giving nontrivial results for 
stochastic efficiency, namely a Mandal-Jarzynski model with 3 states, cf. Fig. 9a. 

It is clear that the amount of delivered work after one entry (trit) of the tape is 
given by: 

w = SE{f-i), (37) 


where i and / are the initial and final state of the system, respectively. The 
corresponding information entropy production in the tape can be defined in multiple 
ways, but if we assume that the dynamics is the same for the reversed trajectory, it is 
given by [8, 18]: 

A stape = k B (In (p/,i) - In (pij)) , (38) 


with i and / the initial and final state of the trit and pj,k the probability for entry state 
k on the tape. The average information entropy change upon processing one trit is thus 
given by: 


(A Stape) =k B ^2 ( Pij - Pfj ) In (pij) . (39) 

j 

Here: 

PFj = e 1 + e «BT +e W j , (40) 

is the probability that the trit leaves the system in state j (j = 0,1, 2). The information 
entropy written in this way is the sum of the change in Shannon entropy and the entropy 
production of an auxiliary Mandal-Jarzynski system if it were to bring the tape back to 
its original distribution [9]. 

The efficiency is defined as the ratio of the delivered work and the amount of 
information consumed. In particular the macroscopic efficiency is given by: 

- = M 

T (As tape ) 

_ <5E (2 (jjp, 2 — Pi, 2 ) + (pf, 1 _ Pi, 1 )) 

k B T ((pi ,0 — Pf,o ) In (Pi, 0 ) + (Pi, 1 — Pf, 1 ) In (pi,i) + (pi ,2 — Pf,z) In (pi, 2 )) 

(41) 

A colour-coded plot is given in Fig. 10, for SE/ksT = 1. Reversible efficiency (fj = 1) 
can be reached in the limit where pij — Pfj, j — 1,2. 

The joint probability distribution Pi n u,N (ni,n 2 ) of n\ and n 2 incoming trits with 
value 1 and 2 respectively on a total of N trits is given by: 


N\ 


(N — n\ — n2)\rti\n2\ 


PiaPiU 1 Pi ,2 


\ N — 711 —712 


(42) 


Plnit,N(nii n 2 ) 
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Pi.i 


Figure 10: Macroscopic efficiency of the Mandal-Jarzynski model in the engine regime, 
with 5E = ksT. The green dot corresponds to pj t i = 0.2, p/ 2 = 0.1 and fj = 0.68. The 
stochastic efficiency for the latter case is represented in Fig. 8 a 




(a) Probability distribution of the (b) Approach of —In P t (y)/t for in¬ 
efficiency. creasing values of t/to (blue, red and 

green curve) to the large deviation func¬ 
tion J(y) (black curve). The yellow 
curve is obtained by extrapolating the 
finite time results. 

Figure 11: Efficiency fluctuations of the Mandal-Jarzynski model, for E = kpT ,pj\ = 
0.2 and p / ;2 = 0.1 . The macroscopic efficiency is given by fj = 0.68 


and the joint probability distribution of n\ and n 2 outgoing trits with value 1 and 2 
respectively on a total of N trits is given by: 

N\ 

PFinal,N{ni, nf) = jtz -TT- j - jP'fiPf^^ ~ PF,1 ~ PFp) " 2 , 

{N — n\ — n 2 )!ni!n 2 ! 

(43) 

which is independent of the distribution of the incoming trits. Using these distributions, 
numerical simulations can be performed, to evaluate the probability distribution of the 
stochastic efficiency rj = w/(TAs tape ), cf. Fig. 11 a. The minimum around 77 = 0 is 
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analogous to the minimum in the effusion model, and disappears in the large time limit. 
Furthermore, one observes the weak local minimum at reversible efficiency 1 with, to 
its right, a more clearly visible maximum. 

The large deviation functions of x\ = rij^/N and X 2 = np 2 /N incoming trits with 
value 1 and 2 respectively, is found from equation (42): 

linit{% 1 , x 2 ) = xi ln(aq) + x 2 ln(x 2 ) + (1 - x 1 - x 2 ) ln(l - x x ~ x 2 ) 

- X\ ln(p 7) i) - x 2 In (p I>2 ) - (1 - x 1 - x 2 ) ln(l - p IA - p /(2 ), 

(44) 

and a completely analogous expression for the large deviation function jFinaiiyi? IJ2) of 
jji = np,i/N and y 2 = npp/N, with pjj- replaced by pF.k- From this, the large deviation 
function of the efficiency can be calculated: 

J(rj)= min (q/nit^l, X 2 ) + iFinaliVl, m)) (45) 

Xl,X2,yi,V2 

where rzq, x 2 , y\ and y 2 are constrained to reproduce the efficiency 17 under consideration, 
i.e., they are related by: 

= SE _ (2 ( 1/2 ~ £ 2 ) + Vi ~ xf) _ 

k B T (x 0 - y 0 ) In (p /)0 ) + (aq - Vi) In (p/,i) + (x 2 - y 2 ) In (p /i2 ) ’ 

This minimisation can be done numerically, cf. Fig. lib. In spite of the noise, which is 
mainly due to the discreteness of variables, the extrapolation seems to work quite well. 
Also, the telltale maximum in the large deviation function close to reversible efficiency 
is again reproduced. 

6. Discussion 

The concept of efficiency plays a crucial role in thermodynamics, especially when 
the efficiency is defined in such a way that is leads to universal system-independent 
statements, such as the one concerning Carnot efficiency. With the advent of stochastic 
thermodynamics, it is natural to revisit such questions for stochastic efficiency. Universal 
statements appear to be possible for the large deviation function characterising the 
asymptotic time regime. In particular, long-time realisations with reversible efficiency 
are exponentially least probable. One purpose of this paper has been to verify and 
document the salient features of the time-asymptotic stochastic efficiency in five different 
settings, namely driven Brownian motion, effusion with a thermo-chemical and thermo¬ 
velocity gradient, a quantum dot and a model for information to work conversion. 

In addition, we provide the analysis for finite time including the approach to and 
extrapolation into the asymptotic time regime. A revealing feature of our analysis is that 
the large deviation properties can be obtained quite consistently by extrapolation from 
rather short finite time results. The other encouraging message is that one can apply the 
analysis to a wide variety of completely different implementations, some of which may 
be easier to realise. Both observations imply that an experimental verification should 
not pose a real problem. In particular, in view of existing experiments on the issue 
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[21, 23, 50, 51, 52], the experimental implementation for the stochastic efficiency of an 
information to work engine should be relatively straigthforward. 
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Appendix A. Large deviation functions from finite time results 

The central object of interest in stochastic efficiency is its large deviation function. 
Unfortunately, the evaluation of exponentially unlikely events is obviously very difficult. 
For this reason, we propose a simple and apparently robust method to deduce the 
large deviation function J{rj) from finite time results. As input we give the efficiency 
probability distributions for three finite times t \, t 2 and f 3 . We propose the following 
ansatz for P t (r/) : 

PmtAv) = Mv)t~ B(r,) e- tJM (A.l) 

Here, A(r t /), Bfq) and Jfq) are three fitting parameters for each r/. J{rj) shall be our 
estimate for the large deviation function. The fitting parameters can be obtained from 
the known values — In (P ti (rj)), with i = 1, 2, 3, since: 


— In (Ptiiv)) /t\ 


1 l/ti In (ti) fti 


J(v) 

-In ( p tM)/h 


1 1/0 In (0)/0 


-In (A(r))) 

_ -ln(P t3 (?7))/t3 _ 


1 I/O In (0) /h 


B(rj) 


Inverting this matrix equality leads to an estimate of the large deviation function. 

Appendix B. Efficiency calculations from event probability distributions 

Consider a model consisting of k tight-coupled processes (e.g. k channels for particle 
transport), where the z-th process undergoes n, events, i = 1.... k. Furthermore, 
the total delivered amount of work and heat can be written as W(ni, ..., nk) and 
Q(ni ,..., n k ). Once the event probability distribution is known, the probability 
distribution of r] can be written as: 

p t(v) = Y (B.l) 

for discrete variables and 

Pt(v) = Jdni..dn k P t (n u ..., n k )8 ^ ) > ( B - 2 ) 

for continuous variables, where P t (ni,..., n k ) is the probability that at time t, for process 
i, n t events have occured. Using the corresponding event large deviation function 
7 (ni,..., rifc) = — lim^oo (1/t In (Pf(ni,..., Ufc))), the efficiency large deviation function 
J{rj) can then be calculated via the contraction principle: 

J(v)= min 7(ni,...,n fc ), 

Til T"i'P'k 


(B.3) 
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where ni, ...,nk are conditioned to: 

-feS 

This provides a good scheme for numerical calculations. 

We proceed to show that there are, in this case, only two extrema of the efficiency 
large deviation function: the macroscopic efficiency and the reversible efficiency. Using 
Lagrange multipliers, we have: 

J(rl) = E*tT Xli ... tXki xL(x 1 ,..,x k ,\), (B.5) 

with 


L(x i,...,xjfe,A) = ...,x k ) + \ (rjQ(xi, ...,x k ) -W(x 1 ,...,x k )). (B.6) 

Therefore, Xi ,..., x k , A are constrained to: 

r-, nri 1 T Al(nr-, nri i \ 

(B.7) 


dn(xi, ...,x k ) _ _ x f dQ(x 1 , ...,x k ) W(x 1 , ...,x k ) 


dxi 


dxi 


dxi 


for all i. As L(x i,..., x k , A) is an extremum of A and aq, i = 1,..., k, we have 

dh(:ri, ...,x fc , A) daq dTQri,..., x k , A) dA dL(xi, ...,x k ,X) 

drj ^ Q x . Qq Q\ Ijq Qq 

_ dL(x i, ...,x k ,X) 
dr/ 

= XQ(xi, ...,x k ). (B.8) 

Note that Q(xi ,..., x k ) = 0 corresponds to reversible efficiency due to the fluctuation 
theorem, and that A = 0 corresponds to d^(x\, ...,x k )/dxi = 0 (using equation B.7) 
which is equivalent with macroscopic efficiency (as we assume convex event large 
deviation functions). As one of these two equalities has to be fulfilled to be in an 
extremum of J(rj), we conclude that these are the only two extremums of the large 
deviation function 

We finally note that the expression for J{rf) can be further simplified if the processes 
are independent. We shall illustrate this for systems consisting of two independent 
processes (with rq and n 2 events respectively). The extension to more independent 
processes is straightforward. The amount of delivered work per event in the itli process 
is written as Wi and the amount of extracted heat is ^, i = 1,2. The efficiency is then 
given by: 

wini + w 2 n 2 


V = 


qi ni + q2n 2 

qini 


r]i + n 2 


Q 2 n 2 


qini + q2n 2 

= ar/i + (1 - a)rj 2 , 


q\n\ + q 2 n 2 


V2 


with 


a = 


qini 


rj-rj 2 


q x n x + q 2 n 2 r/i - r\ 2 


(B.9) 


(B.10) 
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and r]i = Wi/qi, i — 1,2. Therefore, the probability distribution of // at time t is given 
by: 



Pt(v) =11 dwdqP t (w,q)5 ( r) - — 


w 



dn 2 dn 2 "'i''/- - w 2 q 2 \ P t (ni)P t (n 2 )S ( // - m?, - (1 - a)r] 2 ). 

(B.ll) 

From large deviation theory, we can now write the efficiency large deviation function 
J(rj) in terms of terms of the large deviation functions of the event numbers <p, (n,) and 
^ 2 (^ 2 ): 


T/ N , /cuc\ /(I - a)x 

/(,,) = mm ( Vi 


(B.12) 


Appendix C. Event large deviation function of quantum dot-like models 

We present for completeness the large deviation function for the number of particles 
travelling through one energy level of a quantum dot (with energy E and coupling 
constant T = 1) . We only quote the final result, as similar calculations can be found in 
the literature, see e.g. [53, 54, 55, 56]. We recall that the fluxes are given by equations 
(16) and (17). The cumulant generating function is given by: 

/( 7 ) = k + + g(y)s, (C.l) 


with 


k = — 


(kf + k L + k R + k R ^j 


= Jk\k+k L k R - 7o = In ( khkjt 


9 7 ~ V IV \ jL+jL- 

z \ k r k l 

p{ 7) = e 7o/2-7 ; r = k 2 — k:fk R - k+kf; g(q) = p(q) + p(7) _1 - (C.2) 

The large deviation function reads: 


with: 


h(n) = —f( / y(n)) - n^(n), 

2 n / . 

q(n ) = —V r + 2n 2 + 2Vs 2 + rn 2 + n 4 , 

7 c? » 

, , 7o + 4 + q(n) 

7W = y - 


(C.3) 

(C.4) 

(C.5) 
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